Quantum-dot states and optical excitations in edge-modulated graphene nanoribbons 
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We investigate from first principles the electronic and optical properties of edge-modulated armchair 
graphene nanoribbons, including both quasi-particle corrections and excitonic effects. Exploiting the 
oscillating behavior of the ribbon energy gap, we show that minimal width modulations are sufficient 
to obtain confinement of both electrons and holes, thus forming optically active quantum dots with 
unique properties, such as cohexistence of dot-like and extended excitations and fine tunability of 
optical spectra, with great potential for optoelectronic applications. 



Graphene nanostructures have recently triggered a 
wealth of studies for their remarkable properties, which 
combine the unique electronic and mechanical features of 
graphene pQ with the semiconducting behavior induced 
by quantum confinement [2]. Moreover, depending on 
the details of the atomic structure, a variety of novel 
width and edge-related phenomena can arise [3] . To fully 
exploit this richness, great efforts have been devoted to 
achieve precise control of the structure through a num- 
ber of different nanofabrication techniques [4 . Extreme 
control of the width -down to three benzene rings- was 
recently demonstrated by chemical routes for armchair- 
edge graphene nanoribbons (A-GNRs) [5 j, whose optical 
gaps [6-9 finally reach the energy window attractive for 
optoelectronic applications. The perspective of engineer- 
ing ID and 0D quantum confinement of charge carriers, 
so far dominated by unintentional disorder effects [10] , is 
thus becoming realistic [TT] . 

Among the different routes proposed to achieve quan- 
tum confinement in all directions [12], i.e. quantum dots 
(QDs), the most widely used takes advantage of high- 
resolution lithography to carve the full device from 2D 
graphene [13] . This allows to combine atomic-like prop- 
erties -extensively explored in conventional semiconduc- 
tor QDs [14 j- with the advantages born by graphene, 
such as efficient coupling to a graphene-based intercon- 
necting wire or contact and planar geometry compatible 
with available technologies. In such a configuration, both 
dot and barriers are made of graphene by modulating the 
wire width, and the carriers result to be localized in the 
wider region by quantum size effects, as in conventional 
semiconductor nanostructures. 

Here we show that not only quantum size effects, but 
also a novel mechanism -similar to that occurring in con- 
ventional semiconductor hetero] unctions- can give rise to 
confinement of both electrons and holes, where the the 
confining potential landscape is obtained by exploiting 
the peculiar electronic properties of AGNRs. In fact, the 
energy band gap of AGNRs shows three distinct fam- 
ilies depending on the ribbon width, namely N = 3p, 
N = 3p + 1, and N = 3p + 2 (N being the number of 
dimeric lines along the width, and p a positive integer) 
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FIG. 1. (Color online) (a) Atomic structure of the modulated 
graphene nanoribbon. The dot and the barrier regions (N = 8 
and N = 7, respectively) differ in width by one atomic row. 
The armchair edges are H-terminated. (b) Quasi particle 
band dispersions for the infinite GNRs N = 7 and N = 8 (left 
panels), and for the periodically repeated QD (right panel), 
whose unit cell is highlighted in (a) by dotted lines. The two 
non dispersive states (thick magenta lines) are localized in 
the dot (N = 8) region of the structure. Arrows indicate the 
lowest optically allowed transitions (see Fig. [2]). 



[3]: Within each family, the energy gap decreases with 
increasing ribbon width, as expected, but minimal width 
modulations, down to one atomic row, are sufficient to 
induce a large variation in the band gap. For example, 
the conduction and the valence band offsets between the 
N = 7 and the N — 8 ribbon are both as large as 1.4 eV, 
as will be discussed below (see Fig. [I]). 

Joining such ribbons is therefore expected to produce 
an all-graphene system with type-I band modulation sim- 
ilar to a semiconductor interface between different ma- 
terials, while at the same time retaining graphene ir- 
conjugation throughout. In this paper, we provide a real- 
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istic description of QD states and optical excitations for 
such edge-modulated AGNRs by means of state-of-the- 
art ab-initio approaches. Many-body effects, which are 
known to dominate electronic and optical properties in 
low-dimensional systems [7J H5HT7] , are also included via 
the GW and Bethe-Salpeter (BS) schemes [18]. These 
solid-state techniques, though considerably heavy in our 
case, are crucial to study both dot-like and extended fea- 
tures on equal footing, thus allowing us to obtain quan- 
titative predictions of the optical spectra. Our whole 
results demonstrate that edge-modulated AGNRs offer 
a novel mechanism for the creation of optically active 
carbon-based QDs with prominent and tunable exciton 
localization features, which make them suitable for a va- 
riety of applications, ranging from single-photon emission 
to optically-driven quantum information. 

As displayed in Fig. [I] we build our prototype 
graphene-based QD by considering a superlattice ob- 
tained from the periodically repeated junction of two AG- 
NRs belonging to two different families, N = 7 (barrier) 
and N = 8 (dot), so as to maximize the energy gap dif- 
ference and minimize the width variation. Once the two 
constituents of the superlattice are defined, the length of 
both barrier and dot region will then determine the depth 
and the number of quantum-dot states according to the 
confinement mechanism described above. Here we chose 
these parameters in a way to guarantee the presence of 
a couple of localized states maintaining the feasibility of 
calculations [19]. For comparison, N = 7 and N = 8 ideal 
GNRs are also studied. In all systems, dangling bonds at 
the edges are saturated with monoatomic hydrogen [20 . 

The systems described above were fully relaxed by per- 
forming density-functional theory (DFT) supercell cal- 
culations within the local density approximation (LDA), 
as implemented in the Quantum ESPRESSO package 
[2TJ [22] . In order to improve the band structure descrip- 
tion obtained at the DFT-LDA level, we then computed 
the quasi-particle corrections to the Kohn-Sham eigenval- 
ues within the GqWq approximation for the self-energy 
operator. In addition, excitonic effects were taken into 
account by solving the Bethe-Salpeter (BS) equation, 
which describes the exciton dynamics in terms of the 
screened quasielectron-quasihole interaction [18]. From 
the solution of the BS equation, the absorption spectra 
were then computed as the imaginary part of the macro- 
scopic dielectric function. The inclusion of the afore- 
mentioned many-body effects was done using the Yambo 
code [22] [23]. 

Figure [TJb) depicts the quasi-particle band structures 
of the studied systems, that is N = 7, N = 8, and the 
modulated GNR superlattice. In the superlattice (right 
panel), the band gap difference between N = 7 and 
N = 8 acts as a confining potential for the N — 8 region, 
giving rise to an effective one-dimensional potential well 
(periodically repeated). This confining potential is suffi- 
cient to localize two states in the well region (see thick 
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FIG. 2. Optical absorption spectra of N = 7, N = 8, and the 
modulated GNR (QD). The solid black line represents the 
spectrum with the inclusion of the electron-hole interaction, 
while the GW-KPA one (i.e. without e-h interaction) is in 
grey. All the spectra are computed for light polarization par- 
allel to the ribbon axis, introducing a Lorentzian broadening 
of 0.08 eV. 



non-dispersive bands in magenta) , which thus behaves as 
a quantum dot for both electrons and holes. As found for 
ideal GNRs [7] , the quasi-particle corrections to the DFT- 
LDA energy gap are particularly large if compared with 
standard semiconductors, due to a much weaker screen- 
ing and the quasi one-dimensional nature of the systems, 
which both concur to enhance the electron-electron in- 
teraction. This effect is evident also for the superlattice 
structure, where the LDA energy gap between the local- 
ized states is increased from 1.0 to 2.8 eV, that is, the 
GqWq energy gap is almost three times the LDA one. 
In addition, both ideal and modulated GNRs show an 
overall stretching of the band structure of about 20%. 

The optical absorption spectra for the superlattice and 
the two subsystems are reported in Fig. [2j as obtained 
both within the random-phase approximation (i.e. GW- 
RPA, grey shaded area) and including electron-hole in- 
teraction (black line). The main transitions giving rise to 
the low-energy peaks are indicated with vertical arrows in 
Fig.0 As can be seen, the inclusion of excitonic effects, 
which dramatically modifies both peak position and ab- 
sorption lineshape, is crucial to give both qualitative and 
quantitative predictions of all the optical spectra. The 
prominent ID van Hove singularities characterizing the 
RPA spectra disappear giving rise to individual excitonic 
states below the onset of the continuum, with exceedingly 
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large binding energy. This produces an overall red-shift 
of the spectrum, as opposite to the GW gap opening, to- 
gether with a change in the relative position of the peaks 
and possibly in their relative intensity. 

In ideal ribbons, the position of the first peak exhibits 
an oscillating behaviour, according to the family classi- 
fication described above [7]. This allows to span an en- 
ergy window of more than 1.5 eV just by changing the 
width of one atomic row. In the case of the superlattice, a 
peak arises between Aj and Ag, which corresponds to the 
interband transition between the states localized in the 
dot region within the RPA picture [Aqd, see Fig. []Jb), 
right panel]. As previously mentioned, the structural pa- 
rameters of the superlattice crucially determine the num- 
ber and energy location of the confined states, that is, 
the relative peak position wrt A 7 and A$. This would 
further improves the flexibility of this class of systems, 
where a fine tuning of the QD spectrum can be enabled 
both by interfacing different GNRs and by changing the 
dot/barrier length within a given superstructure. 

Let us now better focus on the nature of the first peak 
of the edge-modulated GNR. When excitonic effects are 
included, the first optically active excitation still has a 
predominantly localized nature (i. e. 78% of the weight 
is given by transition between localized states confined 
in the N = 8 region). Thus, the markedly different cor- 
rection to the GW-RPA spectra as compared to ideal 
N=8 GNR has to be ascribed to the additional confining 
potential introduced by the edge-modulation, the first 
one being almost twice as large. This change in the 
confinement properties is also apparent in the excitonic 
wavefunction plotted along the ribbon axis (Fig. |3|: in 
presence of the N = 7 barrier, the overall shape of the 
exciton envelope function changes from Gaussian (top 
panel) to step-like (bottom panel), with a significant re- 
duction of the spatial extension and an exponential decay 
in the barrier region. Nonetheless, the exciton preserves 
a Wannier-like character [Fig. [4^ a)], as found for ideal 
AGNRs 0. 

As described above, while the first excitation has 
mainly a dot-like character, it also contains non negligible 
contributions coming from higher energy levels, thus ac- 
quiring a mixed character. This is true also for higher ex- 
citations, which arise from combinations of single-particle 
states with different localization properties. For instance, 
the second peak is made up of several excitonic states, 
almost degenerate in energy, each of them combining con- 
tributions from single particles states localized in the dot, 
in the barrier, or from resonances delocalized over the 
whole system. The mixed character of the excitations, as 
well as the presence of both dot- and bulk-like excitons, 
have to be ascribed to the unique nature of the system: 
a straddling junction between regions with different en- 
ergy gap is here realized using the same material, thus 
preserving ^-conjugation at the junction interface [24] , 
indeed different from what happens in common hetero- 
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FIG. 3. Electron distribution of the lowest excitons for N — 8 
and the edge- modulated GNR (QD), for a fixed hole position 
(here set to ru — 0). The distribution is plotted along the 
ribbon axis (2), after integrating out the other coordinates. 
For the superlattice ru is chosen to be at the center of the dot 
(N = 8) region. 



junctions made of different materials [25] . 

In order to have further insights in the optical activ- 
ity of such systems, we have investigated the presence of 
dark states in the low energy window for the periodically 
repeated QD structure. Figure [4] illustrates the in-plane 
electron distribution of the the first (bright) exciton (a) 
and few higher optically inactive excitons (b-d) of QD, 
plotted at fixed hole position (black dot). The edge- 
modulations do not introduce optically inactive states 
below the first absorption peak, inheriting the behavior 
of ideal ribbons [7] . Several dark states of distinct nature 
appear instead in the energy region between the first and 
the second peak [Fig. Qb-d)]. In addition to the more 
common dipole forbidden states [see e.g. Fig. |4jb)], we 
find excitons which couple single-particle states spatially 
localized in different regions of the superstructure. Op- 
tical inactivity thus results from the small overlap of the 
electron and hole wavefunctions [see e.g. Fig.[4jc), where 
the hole is in the dot, while the electron is localized in the 
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FIG. 4. In-plane electron distribution of the first bright (a) 
and few higher dark (b-d) excitons for the QD superlattice, 
where the black dot indicates the fixed hole. 
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barrier]. A third type of exciton is shown in Fig. [4^d), 
where the hole is in the dot region and the electron is 
localized in its nearest replicas. 

To summarize, we have investigated the electronic and 
optical properties of edge-modulated AGNRs including 
quasi-particle corrections and excitonic effects. Exploit- 
ing the oscillating behavior of their energy gap, we have 
demonstrated that the simple addition (or removal) of 
single dimeric lines along the ribbon width is sufficient 
to obtain contextual confinement of both electrons and 
holes. We show that these nanostructures can act as opti- 
cally active QDs, whose properties are significantly mod- 
ified by many-body effects. Coexistence of dot-like and 
extended excitations, as well as fine tunability of optical 
spectra are unique features which make these systems 
conceptually different from conventional QDs. These re- 
sults offer a tantalizing perspective, above all in light 
of recent production of atomically controlled armchair- 
edged GNRs. 

The authors thank F. Troiani and G. Goldoni for fruit- 
ful discussions. This work was partially funded by "Fon- 
dazione Cassa di Risparmio di Modena" , Italian Ministry 
of Foreign Affairs, and MIUR-FIRB grant ItalNanoNet. 
CPU time was granted by CINECA through INFM-CNR. 



Corresponding author: deborah.prezzi@unimore.it 
^ Present address: Dipartimento di Fisica, Universita di 
Roma "La Sapienza", P.le Aldo Moro 2, 1-00185 Roma, 
Italy 

[1] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009); A. Geim, Science 324, 1530 (2009). 
[2] M. Y. Han, B. Ozyilmaz, Y. Zhang, and P. Kim, 

Phys. Rev. Lett. 98, 206805 (2007); Z. Chen, Y.-M. 

Lin, M. J. Rooks, and P. Avouris, Physica E 40, 228 

(2007). 

[3] K. Nakada, M. Fujita, G. Dresselhaus, and M. S. Dressel- 
haus, Phys. Rev. B 54, 17954 (1996); Y.-W. Son, M. L. 
Cohen, and S. G. Louie, Phys. Rev. Lett. 97, 216803 

(2006) . 

[4] For review, see e.g. D. Wei and Y. Liu, Adv. Mater. 22, 
3225 (2010); S. Roche, Nature Nanotech. 6, 8 (2011). 

[5] J. Cai, P. Ruffieux, R. Jaafar, M. Bieri, T. Braun, 
S. Blankenburg, A. P. Muoth, Matthiasand Seitsonen, 
M. Saleh, X. Feng, K. Muellen, and R. Fasel, Na- 
ture (London) 466, 470 (2010). 

[6] D. Prezzi, D. Varsano, A. Ruini, A. Marini, and E. Moli- 
nari, Phys. Status Solidi B 244, 4124 (2007). 

[7] D. Prezzi, D. Varsano, A. Ruini, A. Marini, and E. Moli- 
nari, Phys. Rev. B 77, 041404(R) (2008). 

[8] L. Yang, M. Cohen, and S. Louie, Nano Lett. 7, 3112 

(2007) . 

[9] C. Cocchi, D. Prezzi, A. Ruini, M. J. Caldas, and 



E. Molinari, unpublished (2011). 
[10] C. Stampfer, J. Giittinger, S. Hellmuller, F. Molitor, 
K. Ensslin, and T. Ihn, Phys. Rev. Lett. 102, 056403 

(2009) ; P. Gallagher, K. Todd, and D. Goldhaber- 
Gordon, Phys. Rev. B 81, 115409 (2010); M. Y. Han, 
J. C. Brant, and P. Kim, Phys. Rev. Lett. 104, 056801 

(2010) . 

[11] T. Shimizu, J. Haruyama, D. C. Marcano, D. V. 
Kosinkin, J. M. Tour, K. Hirose, and K. Suenaga, Nature 
Nanotech. 6, 45 (2011). 

[12] See e.g. P. G. Silvestrov and K. B. Efetov, 
Phys. Rev. Lett. 98, 016802 (2007); Z. F. Wang, 
Q. W. Shi, Q. Li, X. Wang, J. G. Hou, H. Zheng, 
Y. Yao, and J. Chen, Appl. Phys. Lett. 91, 053109 

(2007) ; H. Sevincli, M. Topsakal, and S. Ciraci, 
Phys. Rev. B 78, 245402 (2008); B. Borca, S. Barja, 
M. Garnica, D. Sanchez-Portal, V. M. Silkin, E. V. 
Chulkov, C. F. Hermanns, J. J. Hinarejos, A. L. 
Vazquez de Parga, A. Arnau, P. M. Echenique, and 
R. Miranda, Phys. Rev. Lett. 105, 036804 (2010); A. K. 
Singh, E. S. Penev, and B. I. Yakobson, ACS Nano 4, 
3510 (2010). 

[13] L. Ponomarenko, F. Schedin, M. Katsnelson, R. Yang, 
E. Hill, K. Novoselov, and A. Geim, Science 320, 356 

(2008) . 

[14] F. Rossi, ed., Semiconductor Macroatoms (Imperial Col- 
lege Press, London, 2005). 

[15] D. Prezzi and E. Molinari, Phys. Status Solidi A 203, 
3602 (2006). 

[16] A. Ruini, M. J. Caldas, G. Bussi, and E. Molinari, 

Phys. Rev. Lett. 88, 206403 (2002). 
[17] E. Chang, G. Bussi, A. Ruini, and E. Molinari, 

Phys. Rev. Lett. 92, 196401 (2004); J. Maultzsch, 

R. Pomraenke, S. Reich, E. Chang, D. Prezzi, A. Ruini, 

E. Molinari, M. S. Strano, C. Thomsen, and C. Lienau, 

Phys. Rev. B 72, 241402(R) (2005). 
[18] For a recent review on GW-BSE see, e.g., G. Onida, 

L. Reining, and A. Rubio, Rev. Mod. Phys. 74, 601 

(2002). 

[19] The length of the N = 7 (N = 8) segment is about 3 nm 
(1.3 nm). For shorter barrier length (e.g. ~ 1.7 nm) the 
first states in valence and conduction band are no more 
localized; enlarging the dot region (e.g. from 1.3 to 2.1 
nm) makes instead the localized states deeper in energy. 

[20] T. Wassmann, A. P. Seitsonen, A. M. Saitta, M. Lazzeri, 
and F. Mauri, Ph ys. Rev. Lett. 101, 096402 (2008). 

[21] P. Gian nozzi et al., \J. Phys. Condens. Matter. 21, 3 95502 

I (2009)| 

[22J SeeEPAPS Document No. XXXX for calculation details. 

[23] A. Marini, C. Hogan, M. Griming, and D. Varsano, Corn- 
put. Phys. Commun. 180, 1392 (2009). 

[24] Type-II nano juct ions with peculiar optical properties can 
also be obtained from AGNRs by means of edge function- 
alization with appropriate chemical groups [9j [26] . 

[25] See e. g. BN-C nanoribbon heterostructures by Sevincli 
et al. [12] . 

[26] C. Cocchi, A. Ruini, D. Prezzi, M. J. Caldas, and 
E. Molinari, J. Phys. Chem. C 115, 2969 (2011). 



